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Abstract 

We calculate the inclusive single-gluon production cross section in high energy onium- 
onium scattering including pomeron loop effects. The resulting inclusive cross section 
is given by the /cT-factorization formula with one of the unintegrated gluon distribution 
functions depending on the total onium-onium scattering cross section, which includes all 
pomeron loops and has to be found independently. We discuss the limits of applicability 
of our result and argue that they are given by the limits of applicability of pomeron loop 
resummation approach. Since the obtained /cT^-factorization formula is infrared divergent 
we conclude that, in order to consistently calculate the (infrared- finite) gluon production 
cross section in onium-onium scattering, one has to include corrections going beyond the 
pomeron loop approximation. 
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1 Introduction 



Recently there has been a lot of renewed interest in resummation of the so-called "pomeron 
loop" diagrams for high energy scattering [T1I21ISII11I3EIIZIIH]- These diagrams are particularly 
important for the calculation of the cross section for high energy scattering of two quark-onia. 
It has been argued that the growth of the onium-onium total scattering cross section with 
energy due to the BFKL evolution is unitarized (at fixed impact parameter) by the multiple 
hard pomeron exchanges realized by the pomeron loop diagrams |im ITU IT^ IT31 IT^ . Such 
unitarization is qualitatively different from the unitarization of total scattering cross section 
in deep inelastic scattering (DIS) on a nucleus [13 CEl UZl UHl UH 120]: in the onium-onium 
scattering case one considers a collision of two small and dilute objects, whereas in the case of 
DIS one of the scatterers is a large and dense nucleus. 

The contribution of multiple pomeron exchanges to the total onium-onium scattering cross 
section has been studied within the framework of Mueller's dipole model ITU] both analytically 
by Mueller and Patel JTj and Navelet and Peschanski and numerically by Salam in |14j . 
The approach of |TT| 118] IT^ has recently been reformulated in the language of the Color Glass 
Condensate by lancu and Mueller in fP.. The main idea behind the dipole model approach 
to onium-onium scattering is reviewed below in Sect. 2. When viewed in the center-of-mass 
frame, each of the colliding onia develops a dipole wave function. The dipoles in the light cone 
wave functions of the onia interact pairwise by two-gluon exchanges fTJ . This process is 
equivalent to resumming pomeron loop diagrams in the traditional Feynman diagram language 
and is illustrated below in Fig. |21 

Our goal in this paper is to calculate the single gluon inclusive production cross section in 
the onium-onium scattering described in the framework of Mueller's dipole model . This 

is an interesting observable one could study in 7*7* scattering at the future linear colliders. 
The calculation turns out to be analogous to the case of gluon production in DIS considered in 
PT] and is outlined in Sect. 3 below. Using the real- virtual cancellations of some of the final 
state interactions observed in ||22j we arrive at the expression for the inclusive cross section 
of the single gluon production given by Eq. (j23|) . As one can also see from a different form 
of Eq. given by Eq. (pUj) . the resulting gluon production cross section adheres to AGK 
cutting rules and is given by the /cT-factorization expression [21], similar to the results of 
[211 12S1 for DIS and to [27|. /cr-factorization appears here because of the fact that, due to 
the nature of the approximation of [TTl , one of the onia wave functions (the one closer to the 
produced gluon in rapidity) appears to be unsaturated, i.e., has no saturation effects in it and 
is described simply by the linear BFKL evolution. This makes the problem effectively identical 
to the case of DIS or pA collisions. 

/c^-factorization formula for the gluon production cross section is known to diverge as ~ 
in the infrared 28j (here kr is the transverse momentum of the produced gluon). This diver- 
gence is beheved to be curable if one includes saturation effects in both of the onia wave 
functions, which would lead to a break-down of the /cT-factorization expression but would yield 
an infrared-finite production cross section (similar observation appears to hold for gluon pro- 
duction in nucleus- nucleus collisions [^ I3U1 1^ ). For the onium-onium scattering, saturation 
in both of the onia wave functions can be included by resumming all pomeron loop diagrams 
with the loops spanning less than a half of the total rapidity interval. However, as we discuss 
in Sect. 4, when rapidity is high enough for such "small" pomeron loops to become important. 
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other corrections to pomeron loops approach become equally important as well. Some of these 
non-pomeron loop corrections have been discussed in jB] to explain the discrepancy between 
the results of |5j and [2,. In the traditional Feynman diagram language these corrections can 
be associated with higher order terms in the impact factors, zero-rapidity size pomeron loops 
and/or with iterations of NLO and NNLO BFKL kernels in the LO BFKL ladder. Therefore, 
we arrive at a somewhat unsettling conclusion that in order to obtain an infrared-finite gluon 
production cross section in onium-onium scattering in a consistent way one may have to go 
beyond the pomeron loop approximation to high energy scattering. 



2 Onium-Onium Scattering at Intermediately High En- 
ergy 

Let us consider high energy scattering of two quarkonia in the approximation put forward in 
IT^ n]. The notations we will use are explained in Fig. ^ The onia have transverse sizes 
2Loi = aio ~" ^1 Xqiii = XQt — Xi, with XQ,Xi,XjQi,Xit the transverse coordinates of quarks and 
anti-quarks in both onia. The impact parameter of the collision is B_ and the rapidity interval 
between the onia is Y. 
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Figure 1: Onium-onium scattering as considered in the text. 

The dipole model description of high energy onium-onium scattering fUElHI is illustrated 
in Fig. 121^.. In the center of mass frame, each of the onia develops a system of color dipoles in 
its light-cone wave function before the collisions. At the time of the collision (denoted by dash- 
dotted line in Fig. |2K) the dipoles in both onia interact with each other pairwise by two-gluon 
exchanges. In Fig. |21A. we show the case when three pairs of dipoles interact with each other. 
This corresponds to three-pomeron exchange in the traditional Feynman diagram language, 
which is illustrated in Fig. |2l3. In the dipole model [TT| IT^. the BFKL evolution [9^, which is 
included in the ladder diagrams in Fig. EB, comes in through the dipole evolution in the light 
cone wave functions of the onia, which leads to generation of the interacting dipoles. 

The dipole wave function of an onium is described by the generating functional for dipoles 
Z{xQi,b,Y,u), which obeys the following evolution equation [TIH IT2] 
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Figure 2: Three-pomeron exchange diagram for onium-onium scattering in (A) dipole model 
and in (B) the Feynman diagram language. The double lines in (A) denote gluons in the 
large-iVc limit. The ovals in (B) denote triple pomeron vertices. (Color on-line.) 
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and m(xoi,6) a function describing a dipole of transverse size at impact parameter b [T2j. 
Using Z{xjQi,b,Y,u) one can define the density of k dipoles in the onium wave function by 

m 



nkixoij h.,y Li, ki, . . . :,Lk-,h.k) 



5^Z{x,,,hY,u) 



(3) 



u=l 



k\ 6u{ri,hi) . . . 6u{rk,h 

If normalized appropriately, would have the meaning of inclusive probability of finding k 
dipoles of given sizes r^'s at impact parameters b^s along with any number of extra dipoles. 
Note that rik, as defined by Eq. (jH)), is normalized to give density of k dipoles rather than 
probability. 

Let us denote the forward amplitude of dipole-dipole scattering mediated by a two-gluon 
exchange as 



d{r,z,b) = (^Ijf In' 



|6+ (l/2)r + (l/2)r' 


\b- (l/2)r - (l/2)r'| 


|6+ (l/2)r - (l/2)r' 


|6- (l/2)r + (l/2)r'| 



(4) 
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Here r and r' are transverse separations of the two dipoles and b is their impact parameter. 

With the aid of the definition from Eq. (0]) we can write the forward scattering amphtude 
for onium-onium scattering, which includes small-x evolution, as 

OO „ 

-^(^01,^0'!',^, Y) = + / (frid\ . . . d^rkd% d\[d\ . . . d'^r[d% 

k=i •' 

X (-l)'=rf(ri,r;,6i-6;) . . . rf(r„r'„ 6, - 6',), (5) 

where we chose the impact parameters of dipoles 01 and O'l' to be B_ and zero correspondingly. 
The notations of Eq. ^ are explained in Fig. CI 
Using Eq. (jS)) in Eq. (0) we re-write it as 

— D(xoi,Xo'i', 5, ^) = — 1 + E y d^rid^hi . . . (frkd%k cPr[d^h\ . . . d'^r'j.d^b'j. 



k=l 



X 



S''Z{x^„B,Y/2,u) 



6^Z{x^,,,,0, -Y/2,v) 



v=l 



X i-lfd{r„r!„b,-I/,) ... d(r„r'„6,-6',), (6) 

where we relabeled functions u by f for the dipole O'l'. If one could solve Eq. ([Q) for an arbitrary 
function u, then the solution could be used in Eq. © to construct the forward amplitude of 
onium-onium scattering. Unfortunately no general analytical solution of Eq. for an arbitrary 
function u is known. 

Eq. (jni) can be recast in a more compact form: 



D{xQi,XQti,,B_, Y) 



exp 



J d^r d% d^r' d%' d{r, r', b - U) 
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(5m (r, 6) (5f (r', 6') 



X Z(xoi,5,r/2,n)Z(xo,i,,0,-r/2,t;) 



(7) 



u=l,v=l 



This is our final expression for the onium-onium forward scattering amplitude here. One could 
also try to rewrite Eq. ^ in terms of a functional integral over functions u and v: however, such 
integral representation appears to require non-Gaussian weight functional, which is probably 
related to the non-Gaussian noise discussed in |2] within the context of adding pomeron loop 
corrections to JIMWLK evolution equations [ini ED] • 

The question of whether -D(xoi, Xq/i/, S, F) gives a unitary total onium-onium scattering 
cross section has been addressed numerically by Salam in It was demonstrated by Monte 
Carlo simulations in {T^ that -D(xoi, ajg'i', S, y) is unitary at a fixed impact parameter b. 
To simplify the analysis of -D(xoi,Xo/i/, 5, F) in Eq. it would be very useful to find an 
integral /differential equation governing the high energy evolution of D{xx)i^Xq,ii^B_^Y)^ similar 
to the case of DIS considered in [TTl HH] . However, despite many attempts, no such equation 
for -D(xoi, Xg/i/, 5, y) has been found at this time. 



4 



Figure 3: Onium-onium scattering amplitude D{Y) from the toy model result in Eq. (|THl 
plotted as a function of agY for d = 0.01 (thick line). The thin horizontal line denotes the 
black disk limit. 



We would illustrate the onset of unitarity in D from Eq. ((7j) by using a simple toy model, 
originally proposed by Mueller in Suppressing the transverse coordinate dependence we 
write the toy model analogue of Eq. (|T)) as |T2] 



dZ 
dY 



asZ'-asZ (8) 
with the initial condition given by 

ZiY = 0,u) = u. (9) 
The solution of the toy evolution equation (jH)) satisfying the initial condition Q is ^T^l 

^(^'«) = ^n"" , .V (10) 
u + [1 — u) e"'' ^ 

The toy model analogue of Eq. is 
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DiY) 



exp 



~d- 



6u 6v 



Z{Y/2,u)Z{Y/2,v) 



(11) 



u=l, v=l 



where d is some positive small number, < c? < 1, and we put Y/2 instead of —Y/2 in the 
argument of Z{Y/2, v) to underline the fact that it is the absolute value of the rapidity interval, 
and not the rapidity itself, which matters there. Substituting Z from Eq. (fTUj) into Eq. (fTTj). 
expanding the exponential and performing the differentiations with respect to u and v yields 



DiY) 



J2 n!(-l)"c;"e°=^(l 



(12) 



n=l 



The series in Eq. fjl2|) is divergent: however, since c? > 0, it is Borel-summable. At large 
rapidities we can approximate 1 — ~ — e"" rewriting Eq. ()12|) as 



D{Y) 



n m jn cts Y 



J2 n\ (-l)"rf"e 



(13) 



n=l 
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which after Borel resummation yields 



- 1-3;^ -P (3^7^)^(0,3^). (14) 

The amphtude D{Y) from Eq. (jl4|) is shown in Fig. |21as a function of as Y. One can see that at 
very large Y the black disk limit {D = 1) is reached and is never exceeded, i.e., the amplitude 
is unitary. 

It is interesting to note that if one does not make the 1 — e"" ~ — e"'' approximation 
and Borel-resums the full series in Eq. fll2|) . the resulting amplitude would violate unitarity: as 
one increases rapidity it would grow, becoming greater than 1 at some large rapidity Y, and, 
after reaching a maximum it would slowly decrease, approaching 1 from above at asymptotically 
high rapidities. Indeed the pomeron loop approximation of jE] |T] consistently resums only the 
leading terms at high energy, and has no control over the subleading terms, such as 1 that 
we neglected compared to e""^/^. Therefore, such pathological behavior of D(Y) is probably 
due to the fact that the toy model at hand, along with the approximation of (TTl [T] , do not 
consistently resum the subleading terms. However, it may also indicate the importance of 
subleading corrections for onium-onium scattering at very high rapidities. After all, the same 
toy model would always give a unitary cross section for DIS, where the dipole-nucleus forward 
scattering amplitude is given by = 1 — Z ^TinHl- We will return to this discussion below in 
Sect. 4. 

3 Inclusive Gluon Production 

To calculate single gluon production cross section in the approximation to onium-onium scat- 
tering developed in ^HQ] we start with the lowest order contribution. One of the corresponding 
diagrams is shown in Fig. EJ The produced gluon is denoted by the cross. The disconnected 
s-channel (horizontal) gluon line in Fig. E] indicates that the gluon can be emitted off either 
the quark or the anti-quark lines in the onium 01. Disconnected t-channel (vertical) gluon line 
denotes the sum of the gluon interactions with the emitted gluon, quark and anti-quark in 
dipole 01 at one end, and with the quark and the anti-quark in dipole O'l' at the other end. 




Figure 4: Lowest order gluon production in onium-onium scattering. 
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The dipole model approach to onium-onium scattering fTT] was developed in the Coulomb 
gauge, which is equivalent to the light cone y4+ = gauge for the wave function of the onium 
moving in the light cone "+" direction and to A_ = gauge for the wave function of the onium 
moving in the light cone "— " direction. In analyzing the lowest order diagram from Fig. |3]we 
will work in the = light cone gauge, where the 01 onium is taken to be going in the "+" 
direction and the O'l' onium is taken to be going in the "— " direction. 

The lowest order gluon production cross section is well-known inSEBllS^- Following 
we can write it as (for a review of similar derivations see [3^) 

do'Lo{XQi,XQ,i,) ^ asCp 1 f ,2 t2 „-ifc-(zi-22) / -I ^1 ~ ^2 ~ 

d^kdyd^B 7r2 (27r)2 J ' \ 

^ dcq (^^1 " 2Ij)2Io'i'5 2^"-*^ ^ ^'^^ (^—2 ~ 2Ij?2Io'i'5 '^—'^ ^ 2ii)^ 



(15) 



where 



^^^^-'-'-^ = ^/^-^'-'^^ = i |6+(l/2)r-(l/2);||fe-(l/2)r + (l/2)rl j ^'^^ 

is the forward amplitude for scattering of an adjoint (gluon) dipole of size r on a fundamental 
(quark) dipole of size r' separated by the impact parameter b and mediated by a two-gluon 
exchange. Again we assume that the impact parameter of the onium O'l' is zero. In Eq. ()15|) . 
k and y are transverse momentum and rapidity of the produced gluon, while z_i and z_2 are its 
transverse coordinates, which are different on both sides of the cut [35J |^ 123 IHE] ■ 

The task now is to include the effects of small-x evolution taking place in both of the 
onia's wave functions into Eq. (jl5p . In our approach we will be closely following the case of 
gluon production in DIS considered previously in |2lj. Just like in Fig. |21 we will consider 
onium-onium scattering in the center-of-mass frame, and, like in |TT] we will be working in 
Coulomb gauge. We will consider the case in which the produced gluon has rapidity y such 
that < y < Y/2, i.e., the gluon is emitted closer to onium 01 than to O'l' in the rapidity 
interval. The rapidity-dependent gluon emissions of Fig. El generating small-x evolution can 
be divided into three categories: (i) emission of "harder" gluons, with rapidity y such that 
y < y < Y/2; (ii) emission of "softer" gluons in the same onium's (01) wave function, i.e., 
emissions in the rapidity interval < y < y; (iii) emissions in the wave function of onium O'l', 
in the rapidity interval —Y/2 < y < 0. 

The analysis of category (i) and (ii) emissions is identical to that for gluon production in DIS 
carried out in Sect. IIIA and IIIB of [21] correspondingly and heavily relies on the real- virtual 
cancellations of the final state interactions discussed in [221 • EI] onium develops a dipole 
wave function which later interacts with the target nucleus. If the onium is moving along the 
light cone "+" direction, the interval of the light cone time x+ over which it interacts with the 
nucleus is negligibly short compared to the time it takes for the dipole wave function to develop 
[10[ 13^ This is very similar to onium-onium scattering: the light cone interaction time 
between the two onia's wave functions is negligible compared to the time required to develop 
each of these wave functions. Thus, from the standpoint of one of these wave functions (the 
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one in onium 01), interaction is instantaneous, happening at, say, light cone time x+ = 0. 
Then, emissions of harder gluons [type (i)] may happen either before or after the interaction, 
at < and x+ > correspondingly. They may be real (the gluon is present in x+ = +oo 
final state) or virtual (the gluon is absent in the final state), given by Figs. 2 and 4 in pT] 
correspondingly. Due to real- virtual cancellations of the final state emissions/absorptions 
of such gluons cancel, leaving only initial state radiation. The overall effect of type (i) emissions 
is to develop a dipole cascade before the collision which would lead to creation of the "last" 
dipole in which the measured gluon is produced. In other words, instead of being emitted by 
the quark and anti-quark in the original onium 01, as shown in Fig. |3] leading to Eq. ()15|). the 
measured gluon would not be emitted in some dipole 23, generated by the dipole evolution. 

The dipole number density at rapidity y in the onium 01 is described by the quantity 
ni{xQi,B_,Y/2 — X23, (l/2)(x2 + ^3)) from Eq. ^ with k = 1. Using Eq. ((T)) we see that 
it obeys the following differential equation tlliil2j, which is actually equivalent to the BFKL 
equation jH]: 



dy 2 TT^ 7 0:20 X 



J d^X2 



20 -^21 L 



^1(^02,^+ ^^21,1/; ^23,^) 



+ ^i(^i2>^+ ^^2o,l/;^23)^) - ni{xQ^,B,y]X23,h) (17) 
with the initial condition 

ni(xoi,5,|/ = 0;x23,6) = (5^(xoi - X23) 5^(5 - 6). (18) 

To include the effects of type (i) emissions one has to modify the cross section in Eq. (fT3j) by 

EH ESI 

daLo{xs^i,Xs^'i') f ,2^ ,2,„ / o v/9 M daLo{x23,XQ,^,) 

d^kdyd^B ^ I dx2sdbn,{x,,,B,Y/2-y,X23,b) ^^^^^^25 ^ (19) 

where we are using the notations of j36j; which are slightly different from those of pT] . 

The emissions of softer [type (ii)] gluons can again be analyzed similar to the DIS case of 
|21j . Analogous to how it was done in Sect. IIIB of [21], one can show that the effect of type 
(ii) emissions is to convert dcq in Eq. (|15|) into the gluon dipole-quark dipole (O'l') forward 
scattering amplitude with a fully evolved wave function of the (upper) gluon dipole. However, to 
complete the picture one has to understand what happens to the wave function of the (lower) 
onium O'l', i.e., one has to analyze type (iii) emissions. This is rather straightforward: the 
lower onium O'l' develops a dipole wave function before the collision. In its light cone time 
x_ the collision is, again, instantaneous, happening at x_ = 0, after which it can again emit 
more gluons. However, none of the gluons emitted in this lower onium wave function is tagged 
upon, i.e., none of them is required to be present in the final state. Therefore, the rules for 
developing the lower onium wave functions are the same as for the calculation of total onium- 
onium scattering cross sections: all final state emissions cancel |221; leaving only the initial 
state emissions shown in Fig. |2K and described by the generating functional Z from Eq. 
[We note also that there is no "cross-talk" (s-channel mergers) between emissions in the upper 
and lower onia in the leading logarithmic approximation considered here.] Combining what we 
know now about type (ii) and type (iii) gluon emissions we conclude that they modify dcq in 
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Eq. (fT3j) into the gluon dipole-quark dipole forward scattering amplitude with a fully evolved 
wave functions of both gluon and quark dipoles. To include type (ii) and type (iii) gluon 
emissions into Eq. (jl5|) we thus have to replace 

doqir, r', b) ^ Dcq (l, r', 6, + |) (20) 

in it. The quantity Dcg can be defined by analogy with D from Eq. (|7j) by noting that, in the 
large- A'^c limit, the gluon dipole is equivalent to two quark dipoles, such that the corresponding 
generating functional is just a square of the generating functional for the quark dipole, 

Zg{xo^,B,Y/2,u) = [Z{x,„B,Y/2,u)f, (21) 

so that 



DGq{xQi,Xo,i,,B,Y) = 1 - <^ exp 



- J (fr d% d^r' d%' d{r, r', b - U) 



6u{r, b) 5v{r!,y 

X Z(xoi,5,r/2,M)Z(xoi,5,r/2,M)Z(xoa',Q,-r/2,t;)| . (22) 

/ u=l,v=l 

Combining Eqs. (jl9p and (j2m) we finally write for the inclusive gluon production cross 
section in onium-onium collisions 



da{xQi^XQiii] 
d^kdyd^B 



J d^X2d^X3 ni (^Xoi,5, ^-y;x23, ^(^2+^3)) 



X 



(27r)2 7 kl-Xi\^ k2-Xj 

1 Y\ / 1 Y 

Dcq ( zi - Xj,Xo,i,, -{z^ + Xj),y + —j + Dcq U2 - ^i,^o'i'> ^{z2 + Xi),y+ — 

I, . Y\ ^ / 1. . F 



Daq [Zi- Z2,Xo,y,-{z^+ Z2),y + —j - Dcq [^i - Xj,Xjj,y, -{Xi + Xj),y + - 



(23) 



This is our main result. The above formula is valid for < y < Y/2, and can be easily modified 
to give the production cross section for —Y/2 < y < by switching 

y ^ -y, xoi ^ XQ/r. (24) 

If y = (gluon is produced at mid-rapidity) Eq. (j23|l becomes symmetric, since Dgq should 
contain only the linear BFKL evolution, due to the fact that the approximation to onium-onium 
scattering employed here d does not consistently resum pomeron loops with rapidity extent 
less or equal to Y/2 (see Fig. El below). While formally Eq. with y = sums up some 
pomeron loops of rapidity size Y/2 and smaller in Dcq, as given by Eq. (j221), such loops were not 
consistently summed over in arriving at Eq. (j^^ . Therefore, all but the linear BFKL pomeron 
terms in Dcq should be neglected in Eq. ((221) at ?/ = 0. 

It is interesting to notice that the evolution between the produced gluon and the closest 
onium (onium 01 in Eq. ((221)) is given by the linear BFKL equation. This result is in accordance 
with the AGK cutting rules [221 ■ This is indeed surprising, particularly due to the fact that 



9 



AGK violations have been found recently in inclusive production of qq pairs jSZI and in inclusive 
two-gluon production JHH] in DIS and pA collisions (see pB] for a relevant discussion). However, 
the reason our result (j23|) adheres to AGK rules is the same as for the single inclusive gluon 
production cross section in DIS: due to the nature of the approximation to the onium-onium 
scattering which we employed here following |1H , there is no difference between the wave 
function of the onium closest to the produced gluon in rapidity and the dipole wave function 
of the qq pair in DIS, considered in [21]. Both wave functions only include pomeron splittings. 
This leads to the similarity between Eq. and Eq. (30) in |21jj and to the fact that both 
cross sections adhere to AGK cutting rules. While we can not explain why some inclusive cross 
sections [21] adhere to AGK rules and some other cross sections jSZl EH] violate them, it is 
clear that the reason AGK rules work for onium-onium scattering as well as for DIS is due 
to moderately high energy approximations applied with the non-linear saturation effects being 
included only in one of the scatterers wave functions. 

Similar to the result of j2I], Eq. ()23|1 can be recast in fc-r-factorization form |21]. First we 
note that, as long as impact parameters are concerned, ni{x,B_,Y]r,b) depends only on B_ — b. 
Therefore, 



r / Y 1 



Y 



(25) 



in Eq. depends on X23 = ^2 ~^3) but does not depend on 212 +213. Using this fact we can 
recast Eq. (j^!^ into the following form, similar to how it was done for Eq. (30) in (note 
that now we do not assume that the amplitude Dcq is independent of impact parameter, as we 
did in 



•j 



X 



(Pkdy (2-nf kf 



p j d^z S/l j d'^bDGg(^z,Xo,y,b,y+^^ 



-ik-z 



z — Xij) In 



- ■ Til I Xgi, — y, X23 

\z — XiAiil \ 2 



(26) 



i,j=2 - \\^ 

For derivation of Eq. see Appendix A. Defining the unintegrated gluon distribution func- 
tions in the onia by 

3 



as (27r)= 



X23 (z - In 



1 



ij=2 



I ^ 3y 2 -I I 



X ^1 (^01.1/; ^23) 



and 



4>o'i'ik,y) 



j d''zd%e-''^-^VlDGq {z,x^,,,,b,y) 



as (27r)3 

we can rewrite Eq. (f^Bj) in the familiar fcT^-factorized form |211 123 1^ 

dcr{xj[)i,XjQiii) 2 a., 1 f ,n . / Y 



d^kdy 



(27) 
(28) 

(29) 



We have thus shown that the inclusive gluon production cross section in onium-onium scattering 
analyzed in the approximation of fn\ IT] is given by the fc-r-factorization formula ((23) ■ 
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4 Limits of Applicability of the Pomeron Loop Resum- 
mation 



Eq. (j^^ presents a problem. As long as 0oi and 0o'i' are positive definite, the inclusive cross 
section in Eq. (j29|) diverges as ~ 1/]^ in the infrared (at low kx = \k\)- This is a well-known 
problem of perturbative gluon production [221 ISSl IMj resulting in the infinite total number 
of gluons produced. The divergence is usually ~ at the lowest order (for pp collisions) 

|32[ IHHl El] shown in Fig. |3] It reduces to ~ l/k"^ for pA collisions and DIS [33] , when the 
saturation effects are included in the target nucleus. There is a strong belief in the community 
that including saturation effects both in the target and in the projectile wave functions would 
eliminate the power-law divergence completely. (Elimination of divergence would indeed require 
the fc-r-factorization formula to break down, since no matter what the gluon distributions are it 
would always lead to ~ l/k"^ singularity in the IR. ) The belief is confirmed by numerical [301 
and analytic analyses of gluon production in AA collisions. Therefore, one would expect 
that a complete description of onium-onium scattering should also lead to infrared-regular gluon 
production cross section. The fact that our result does not give such cross section, indicates 
that the approach to onium-onium collisions from jTT| |Tj applied here is incomplete. (Indeed 
the same problem is present in DIS and pA collisions in the approach used in jTTl ESI EI]-) 

To see what is missing in the approach of pTf llj let us perform a few parametric estimates. 
First of all, a single BFKL pomeron exchange over rapidity interval Y is, parametrically, of the 
order of 

a^e^-^-i)^, (30) 

where 

ap-l = — In 2 (31) 

TT 

is the intercept of the BFKL pomeron |9|. The exponent in Eq. (jHUI) comes from ladder rungs 
resummation with the factor of in the traditional Feynman diagram language coming from 
coupling of the ladder to the onia (impact factors). In the dipole model the factor a1 comes from 
the two-gluon exchange interaction between the two dipoles at the ends of dipole evolutions in 
both onia [Eq. (jH)]. 

Pomeron loops become important in onium-onium scattering at rapidity Yu at which the 
correction ()3U|) brought in by an extra ladder becomes of order one, i.e., when 

^2g(ap-i)y^ ~ 1, (32) 

leading to [T^ITT] 



^c/ ~ T 7ln— . (33) 



1 1 1 
7 In — 

ap — 1 

Therefore, for rapidities Y>^Yu pomeron loops become important and have to be resummed. 

The types of pomeron loops included in the approach proposed in [TTJ |T] are shown in 
Fig. EK. There the thick lines denote BFKL ladders (hard pomerons). The approach of |1H IT] 
allows only for pomeron splittings at positive rapidities (above the dashed line) and only for 
pomeron mergers at negative rapidities (below the dashed line). Obviously this approximation 
consistently resums all pomeron loops stretching over rapidity interval greater than F/2, where 
Y is the total rapidity interval. However, some of the "smaller" pomeron loops, namely the 
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Figure 5: The pomeron loop diagrams included (A) and not included (B) in the approach to 
onium-onium scattering of . Thick lines denote BFKL ladders. 



loops stretching over rapidity intervals less than Y jl^ are not included consistently. An example 
of pomeron loops which are no included is shown in Fig. 03, where the pomeron ladder both 
splits and merges at positive and at negative rapidities. One may argue that small pomeron 
loops are always subleading at high energies and can be safely neglected compared to larger 
loops. Still, the small loops of Fig. 03 become order 1 corrections at rapidity defined by 

^2g(ap-i)y^/2 _ (34) 

which gives 

Yv = 2Yu. (35) 

Therefore, we conclude that large pomeron loops of fn\ are a good approximation to high 
energy onium-onium scattering cross section only in the rapidity interval 

Yu <Y < 2Yu. (36) 

But what happens at higher rapidities, for F > 2 Yjy? Would resummation of smaller 
pomeron loops of Fig. |3[3 be sufficient to consistently calculate the total onium-onium scat- 
tering cross section? Our answer is "no". The reason for this conclusion is that other (non- 
pomeron loop) higher order corrections also become important ai Y > 2Yu. For instance, 
NLO corrections to pomeron impact factor lead to a contribution of the order of |39j 

«3e(-^-i)^ (37) 

coming from the whole ladder. One may also speculate that inserting one rung of NLO BFKL 
kernel ^ into LO BFKL ladder jJT] would generate the contribution of the same order as 
shown in Eq. (jTTj) . The corrections from Eq. (|H7j) would become of order 1 at rapidity (3/2) Y[/. 
This would lead to an even smaller upper bound in Eq. which would become (3/2) Y"^/ 

instead of 2 Yu. The resulting applicability window for large pomeron loops is then limited to 

Yu<Y<^Yu. (38) 
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Indeed, since y^/ 3> 1 when <^ 1, the rapidity window in Eq. (j38|) would still be very large 
for asymptotically small coupling. Note that, for the values of the strong coupling constant 
accessible in the current accelerators, NLO BFKL corrections are numerically large and would 
modify our above conclusions. 

At rapidity Y = 2Yu other corrections will become important as well. In performing the 
integration over rapidities of a leading order pomeron loop the lower limit of integration gives 
a loop stretching over zero rapidity interval fSl ^| • Such diagram would still have a single 
pomeron exchange, giving e^"^"^-*^, and the zero size loop, giving due to two pair of dipoles 
interacting, overall contributing a factor of 



■1)Y 



(39) 



A correction like the one shown in Eq. ()39|) may also result from inserting NLO BFKL kernel 
|in] twice in the LO BFKL ladder [H], or from inserting NNLO kernel once in the LO BFKL 
ladder. 

The corrections in Eq. (j39|l become of order 1 when Y = 2Yu. In other words, these terms 
become comparable to "small" pomeron loops stretching over Yjj rapidity interval. This is 
natural since the resummation parameter for these higher order corrections is easily expressed 
in terms of small pomeron loops resummation parameter 



c,4e(«P-i)y 



^2g(ap-l)y/2 



(40) 



so that both become comparable to 1 simultaneously. 

We can conclude that pomeron loop approximation is only valid in the rapidity region 
specified by Eq. flSHjl . To go beyond that rapidity interval one has to include a scope of other 
subleading non-pomeron loop effects listed above. The corrections of Eq. ()39p were recently 
discussed in jH] to account for the apparent discrepancy between the results of j21 and 
indicating that such corrections can not be consistently resummed within the current pomeron 
loop resummation formalisms. 

Now let us return to the question of gluon production in onium-onium scattering. A typical 
diagram contributing to the above result ()23|) is shown in Fig. There, the produced gluons 
are denoted by outgoing double arrows and the BFKL ladders are again denoted by thick lines. 
The graph in Fig. IHK has linear BFKL evolution (single ladder) on one side of the produced 
gluon (in rapidity), with all the nonlinear pomeron loop effects being on the other side of the 
produced gluon, just like we have obtained in deriving Eq. (f^Hj) . The resulting cross section is 
IR divergent as ~ and is valid in the rapidity window given by Eq. (jHHj) . Therefore, the 

^ divergence persists as long as pomeron loop approximation is valid. 

Small pomeron loop corrections of Fig. 03 become important at rapidities Y > 2Yjj. They 
may lead to gluon production diagrams like the one shown in Fig. Such graphs include 
non-linear saturation effects in the wave functions of both onia. However, they do not violate 
fcr-factorization, as could be seen in the graph of Fig. IHfB, which has no "cross-talk" between 
the onia wave functions other then the interaction leading to gluon production. Therefore, 
diagrams like the one shown in Fig. would still lead to ~ l/Zcf^ divergence in the infrared. 

However, one may also imagine small pomeron loop corrections which would violate kx- 
factorization, as shown in Fig. EP. As we discussed in Sect. 3, in the large pomeron loop 
approximation of fn\ such corrections would cancel due to real- virtual cancellations of 
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ABC 

Figure 6: Gluon production diagrams for onium-onium scattering: (A) a typical diagram con- 
tributing to Eq. (j^^ above; (B) contributions of smaller pomeron loops still giving IR-singular 
gluon production cross section; (C) contributions of smaller pomeron loops which may lead to 
IR-regular gluon production cross section. Again thick lines denote BFKL ladder and outgoing 
arrows denote produced gluons. 

making the final result for the inclusive cross section adhere to AGK cutting rules. Once the 
small pomeron loops of Fig. are included in both of the onia wave functions, such real- virtual 
cancellations may stop being valid. This would lead to diagrams in Fig. contributing to 
gluon production cross section, and, possibly, regulating the ~ divergence. However, 

such regularization of the divergence would be due to small pomeron loops, which become 
important at rapidities Y >2 Yu, when the pomeron loop approximation breaks down. Thus it 
may happen that regularization of ~ divergence would have to also include effects going 
beyond the pomeron loop approximation, alongside with the purely pomeron loop diagrams of 
Fig. Ep. 

If, instead of going to higher energies/rapidities, we fix the rapidity Y and decrease kr, 
"higher twist" effects would become progressively more important. Our conclusion above im- 
plies that, while higher twist terms due to large pomeron loops of [TTl [T] do not regularize 
the ~ divergence of the inclusive gluon production cross section, the higher twist terms 
associated with small pomeron loops, not included in the approximation of |1H IT], which may 
regulate the ~ l//c|. singularity, come in at the same time as the higher twist terms associated 
with other (non-pomeron loop) kinds of corrections. 
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Appendix A 



Starting with the first term in the square brackets of Eq. ()23|) (integrated over B) we note that 
Z2 integral can be carried out using 



yielding 



dhe'^-^^ = 2 7Tt^ (Al) 



\—\ —i I ^ 



/ 1 Y\ 

X Dcg Ui - xj , xoa' , - Ui + ^j) , y + — j • (A2) 

Defining z = Zi — Xj and rewriting d'^X2CPxs = d'^X23(fb, where b = Xj + {l/2)z we can 
transform Eq. ()A2|) into 



'y . . I ^ 



X J d%DGq(^z,Xo,r,b,y + —y (A3) 

The second term in the square brackets of Eq. (j23|) gives an identical contribution. 
The third term in the square brackets of Eq. is different. Defining 

— 1 — —1 —i 

and 

i.2 = ^2 ~ 3ij W6 write it as 

Y \ a,C 



j d^X2 d^xs ni fxoi, — - y, X23) / ^ / d^zi d^i2 ^ ■ ^ J2 ("l)* 
J \ 2 J [2 Try J z^ z^ 



1 . . 1 Y 
-{zi + is) + 2 + 2/ + 

Dropping tildes over z^ and rewriting (i^X2 d'^x^ = d'^X23 d% with 



X Dcq (zi - Z2 + Xij,XQ,y, -(^1 + Z2) + -{Xi + Xj), y + — )- (A4) 



b = ^{Xi + Xj) + ^(^1+^2) 



we recast Eq. (IA4|) into 



- [ d^X2z rii (xq^, — - X23) f [ d^zi d^Z2 ^ ■ ^ J2 ("1) 
J \ 2 J {2 7ryn^ J zf zi f-^^ 



Y^ 



X J d%DGq(^Zi- Z2 + Xij,Xo,i,,b,y + . (A5) 
Defining z = z_i — Z2 + x^j and integrating over y = (1/2) (^^ + Z2) using 



j2 y y + ^ ^1 1 



^^2 -^^ = 2^1^^ (A6) 
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with /i some infrared cutoff yields 



(A7) 

Finally, similar to the above, the last term in the square brackets of Eq. ()23|) can be written 



as 



- / d^X23 ni (xoi,^ -1/;X23 



3 



Combining Eqs. jSHl) (twice), and (fA8|l yields 



)^0'1' 



2 7r2 



which can be rewritten as 

da{xjoi,Xjo,i,) ^ 
d^kdy 



Y 



2ik- - -'\ 

\z — XiA 



- In 



p / E (-1) 

- 271 5^{z-x 



(A9) 



d^X23 rii (xoi,^ -i/;x23 



OsCf 1 

2 7r3 F 



i,j=2 



e-'^-^ In 



\ I ^ i j I / 

where is the transverse coordinate (2d) gradient squared. Integrating by parts yields 



(AlO) 



d^kdy " 2 7r3 F 

X E (-i)^+^' / rf^ 

and, remembering that 



d\ 



X23 ni ^01 



Vi I d% Dcq [z,xj^,y,h,y + 



Y 



-ik-z 



VM^Mn- 



y] X23 In 



4 In 



I 3/ ^' ^' I ^(X 



1^1 



1^1 



(All) 



(A12) 



reduces Eq. (jAllj) to Eq. as desired. 
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